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Abstract 

X. 

^ , We analysis the symmetries of the reflection equation for open XY Z model and 

find their solutions case by case. In the general open boundary conditions, the 
Lax pair for open one-dimensional XY Z spin-chain is given. 
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1 Introduction 



In the formalism of the quantum inverse scattering method (QISM), the integrabihty of 
a system with periodic boundary conditions is based on the Yang-Baxter equation|l|, 0. 
Of central importance in the QISM is the 7^-matrix. Relying on the previous work of 
Cherednik [Q, Sklyanin generalized the QISM to system on a finite internal with inde- 
pendent boundary conditions on each end]^. For such a system, the boundary conditions 
are defined by the operators K"^, which satisfy the so-called reflection equations (or called 
the boundary Yang-Baxter equations). The integrabihty conditions for systems with open 
boundary conditions are that both the Yang-Baxter equations and the reflection equations 
must be satisfied. Since Sklyanin's work, a lot of integrable models with open boundaries, 
both in quantum field theory and in statistical mechanics, as the Toda chain[0, nonlinear 
Schrodinger equation[|[, the A2 spin chain with Uq[SU{2)] symmetry]^], supersymmetric 
t-J modell, 0, sine-Gordon modelg], v4„„i vertex models[|], the XYZ modelg, |10|, |lT 
have been studied extensively. 



In this paper, we will present the solutions K'^{u) to the reflection equations of open 
XYZ model. The solutions K^{u) for XYZ model have been studied by de Vega, 



Gonzalez-Ruiz and lnami[p!0|, ^A\ . However, we will persue a new procedure to the prob- 
lem. Firstly, we analysis the symmetries of the component forms of the reflection equations 
for XY Z model and determine the equations to be solved. Then we study the chosen 
equations to get the solutions case by case and their corresponding boundary terms 
in the Hamiltonians. For the case of trK^iQ) = 0, we get a new K{u) that is not given 



in ref. [|T0[ (see (2.18) and (2.25)). Furthermore, we find that several if-matrices may cor- 



respond to the same boundary terms if they differ by a scaling function satisfied some 
conditions. We expect our procedure and discussions to be helpful to the understanding 
of the role played by fC-matrices in the description of a system. In the formalism of QISM, 
the systems described by the above Hamiltonians with boundary terms are integrable. 

However , in the case of periodic boundary condition, there is an equivalent approach. 
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which is called the Lax representation 0, to the proof of the integrability of the system. 
In this formalism, a model is said to be completely integrable if we can find a Lax pair 
such that the Lax equation is equivalent to the equation of motion of the model. Now, 
it is natural that one expect a variant of the Lax representation to systems with open 
boundary conditions. In fact, this has been done, for example, for the one- dimensional 
(ID) Heisenberg open XXZ chain |jl2[, the Hubbard model|jl3[. We will construct the 
Lax pair for the ID open XYZ spin-chain. 



The paper is organized as follows. In section 2, after briefly reviewing the 7?.-matrix for 
the XYZ model, the corresponding Yang-Baxter equation and the reflection equations, 
we introduce the symmetries for the reflection equations and find the i^'-matrices case by 
case. In the second half of this section, we give the corresponding boundary terms for 
each solution and discuss a symmetry of the general Hamiltonian for the open XYZ 
model. In section 3, we will construct the Lax pair for the XYZ model with general 
boundary conditions. We will make some discussions in the last section. 



2 Solutions to the reflection equations for the XYZ 
model 

The XYZ model is defined in terms of the Boltzmann weights given by the elliptic eight- 
vertex solution of the Yang-Baxter equation: 



7^12(M)7^13(M + v)n23iv) = 7^23(^^)7^13(M + v)ni2iu) 



(2.1) 



Here TZu^u) = lZ{u) ® 1,7^23(m) = 1 <S)'R.{u) . The 7^-matrix 7l{u) is given as follows: 

^ sn(M + 1]) ksnrjsiausn^u + ?]) ^ 

snu snr] 

snr] snu 

fcsn?7snMsn('U + rj) sn{u + f]) J 



n{u) 



(2.2) 



where sn (and cn,dn in the formulas below) is the Jacobi elliptic function of modulus 
< k < 1. For the properties of the 7^-matrix (2.2) see refs. |I0|, |Tl| . 



To find boundary terms compatible with integrability, Sklyanin has introduced a pair 
of boundary /^-matrices if^lQ. are the solutions of the following reflection equations: 

7^12(M - v)K{{u)'JZi2{u + v)K^{v) = K^{v)ni2{u + v)K{{u)'}Zi2{u - v) (2.3) 

ni2{-u+v)K+''{u)ni2{-u-v-2r])K^''{v)=K+''{v)nu{-u-v-2r])Kt''{u^^ 

(2.4) 

where Kf" = 1, K2 = 1 K^. Here we should note that equations (2.3) and 
(2.4) only hold for the symmetric 7?.-matrix(i.e. TZ has both P and T symmetry). The 
7^-matrix given in equation (2.2) is symmetric. 

For a solution of the equation (2.3), then 

K+{u) = K-\~u-7]) (2.5) 

gives the solution to the equation (2.4). So it is sufficient to consider the equations (2.3). 
Let K~{u) be of the following form: 



K-(u) 



^ ai{u) a2{u) 



(2.6) 



For the 7^-matrix (2.2), the equation (2.3) in component forms is equivalent to the fol- 
lowing 12 equations: 

die2U2V2 + C1U3V2 - C1U2VS - die2UsVs = 0, 
^261^2^2 + C2U3V2 - C2U2VS - d2eiUsVs = 0, 



(2.7) 



e2UiVi + eiu^vi — eiUiv^ — 62^4^4 = 0, 

Cid2UiVi + C2(iiM4f 1 — C2(ilMlf4 — Cid2U4V4 = 0, 



(2.8) 



C2eiU2fi - Cie2U2Vi - diu-iVi + d2U3Vi - eiUiV2 - 62^4^2 + C2diUiV'i + cid2U4V-i = 0, 

0261^2^4 - Cie2U2V4 - diU^V^ + d2U3V4 - eiM4f2 - e2Mlf2 + C2diU4V3 + Cid2UiV3 = 0, 

(2.9) 



C2eiU3Vi - Cie2U3Vi - diU2Vi + d2U2Vi - eiUiV3 - e2U4V3 + C2diUiV2 + Cid2U4V2 = 0, 
C2eiU3V4 — Cie2U3V4 — diU2V4 + d2U2V4 — eiU4V3 — e2UiV3 + C2diU4V2 + Cid2UiV2 = 0, 

(2.10) 



d2eiU2Vi + C2U3V1 — U1V3 + did2UiV3 + C1C2U4V3 — eie2W4'J^3 - die2U2V4 — C1U3V4 = 0, 

d2eiU2V4 + C2U3V4 — U4V3 + did2U4V3 + C1C2U1V3 — eie2Mi'i;3 - die2U2Vi — C1U3V1 = 0, 

(2.11) 



d2eiU3Vi + C2U2V1 - U1V2 + did2UiV2 + C1C2U4V2 - eie2M4'y2 - die2U3V4 - C1U2V4 
d2eiU3V4 + C2U2V4 - U4V2 + did2U4V2 + C1C2U1V2 - eie2UiV2 - die2U3Vi - C1U2V1 

where Ui — ai{u),Vi — ai{v){i — 1, 2, 3, 4) and 



= 0, 

= 0, 
(2.12) 



sn{u+v+ri) ' 
sn{u+v) 



C2 

do = 



snr? 



sn{u—v+ri) ' 
sn{u—v) 



sii{u+v+r]) ' ^ sn{u—v+Ti) ' 

Cl — ksnr)Sii{u + v) , e2 — ksn'r)sn{u — v) . 



(2.13) 



After analyzing the system of the equations (2.7) -(2.12), we find that there are three 
kinds of symmetries for the solutions of equation (2.3). 

(A) The Scaling Symmetry. Multiplication of the solution K''{u) by an arbitrary function 
f{u) is still a solution of equation (2.3). 

(B) Symmetry of Spectral Parameter. If we take a new spectral parameter u = fiu where 
/X is a complex constant, the new matrix K~{u) is still a solution of (2.3). This symmetry 
is useful when we consider the rational limit of the matrix K~ (u) . 

(C) Symmetry of Interchanging Variables u and v with Different Indices. Under the ex- 
change of U2 ^ U3 and f2 fa, the equations (2.7) are invariant ,while the equations 
(2.9), (2.11) are turned into the equations (2.10) and (2.12), respectively. On the other 
hand, upon the exchange of ui U4 and vi <->■ V4, the equations (2.8) are invariant, while 
for each set of equations (2.9)-(2.12), one equation is changed into another. 
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In view of the relation among Ci,di,ei{i = 1,2) 

die2 d2ei 



(2.14) 



Cl C2 

two equations in (2.7) and (2.8) are equivalent, respectively. Because of this fact and the 
above symmetries, we only need choose four equations from equations (2.7)-(2.12) in order 
to solve the reflection equation (2.3). In the following, we will take the first equation of 
(2.7), (2.8), (2.9) and (2.11), respectively. Of the four chosen equations, three equations 
are necessary to determine the K~ due to the symmetry (A). The fourth equation may 
be as the consistency condition for the solutions of the three equations. 

The solutions can be classified according to whether ai and 02 equal to zero or not. If 
assuming 02 7^ 0, we differentiate both sides of the first equation of (2.7) with respect to 
the spectral variable v and set u = v, then we get the following equation 

^{^) = -ksni2u)[l - (2.15) 

Similarly, assuming oi 7^ O,one also has 

i_(!fl) = (2.16) 

if we differentiate both sides of the first equation of (2.8) with respect to v and set u = v. 
Now we solve the four equations chosen above case by case. The procedure we adopt 



below is the same as that in refs.[|Tn, O]. We give the main results in the following. 



Case (a): ai = 0. Using the first equation of (2.8), we have 04 = 0. The nontrivial 
cases correspond to 02 7^ and 03 7^ which come from the first equation of (2.7). It has 
the following two subcases. 

Subcase (al): (^)^ = 1- From equation (2.15) and the symmetry (A), we obtain 




K-iu) = (2.17) 
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Subcase (a2): (^) 7^ 1- Integrating equation (2.15) and taking account of the symmetry 
(A), we have 



K-(u) 



\ 



A_(l — ksn^u) — 1 — ksn^u 
here A_ is an arbitrary constant. 



(2.18) 



Case (h): 02 = 0. Prom the first equation of (2.7), we have — 0, so the nontrivial 
cases are those of ai 7^ and 04 7^ from the first equation of (2. 8). Here we also have 

two subcases. 



Subcase (hi): (^)^ = 1- From equation (2.16) and the symmetry (A), we get 



K-{u) 



1 
±1 



(2.19) 



Subcase (h2): (^)^ 7^ 1- Considering the symmetry (A) and integrating equation (2.16), 



we have 



K-{u) = 



sn^_ 



sn(^_ +u) ^ 

sn({_ - u) J 



(2.20) 



here ^_ is an arbitrary constant. 



Case (c): Oi 7^ 0, 02 7^ 0. Using the first equation of (2.7) and (2.8), respectively, we 
know that 03 7^ 0, 04 7^ 0. For this case, we may have following four subcases. 
Subcase (cl): = 1, (^)^ = 1, i.e. M3 = £2^2,^4 = £1^1, — el — 1. Using the first 

equation of (2.9), we get two K~ solutions: 



1 + e2ksn^u tt_sn-u 
£2Q!-snM 1 + £2^sn^ti 



{for £1 = 1) 



(2.21) 



{for £1 = -1) 



(2.22) 



1 + £:2^sn^-u /9_cn-udn-u 
y e2/9-cntidnti —1 — S2ksiL^u 
where q;_,/3_ are two arbitrary constants. 

Subcase (c2): {f-f = 1, {^f ^ l,Then 1^3 = e2U2,el = 1. Solving the first equation of 



(2.9), we have the following solution 
1 I 

K-iu 



sn^_ 



a_snticnMdnM 



V 



sn(e-+M)(l + £2fen2M) i-fc2,,2g_,,2. 



(2.23) 



1— fc^sn^J— sn^n 

here ^_ and «_ are two arbitrary constants. 

Subcase (c3): (j^)^ 7^ 1, (^)^ = 1- Taking M4 = and solving the first equation of 
(2.9), we get two solutions, 

/i_sn(2-u)[(l — A;sn^-u)A_ + 1 + fcsn^-u] 



cnwdn-u 

/i_sn(2'u)[(l — ksn^u)X- — 1 — ksn^u\ 
{for El = 1) 



cnudnu 



(2.24) 



snu 

/i_sn(2'u)[(l — A;sn^M)A_ 
(for El = -1) 



/x_sn(2-u)[(l — A;sn^-u)A_ + 1 + ksir^u] 
ksn^u] —snu 

(2.25) 



here A_ are two arbitrary constants. 

Subcase (c4): (^)^ 7^ 1; (^)^ 7^ l.This is the most general case. Now we first have to 
integrate equations (2.15) and (2.16), then insert the results into the first equation of 
(2.9). Upon the use of the symmetry (A), we find that K~ has the following form: 



/ 



K-(u) 



sn^_ 



sn(^_ + u) 



/i_sn(2u) 



(1— fcsn'^u)A_ — 1— fcsn^u 
1— fc-^sn^^-sn^u 



/n \ (1— A;sn^'ii)A_+l+fcsn^« 
U__Sn(2'U)-^^ , ,2 2/: — 2 

sn(^_ - u) 



(2.26) 



here A_ are arbitrary constants. 



Note that the subcases (al),(bl) and (cl) are discussed in ref.[|^, while the subcase 



(c4) is studied in ref.pT 



Now we consider the boundary terms corresponding to the i^-matrices given in equa- 
tions (2.17)-(2.26). For the case of trK^{0) ^ 0, the Hamiltonians with boundary terms 
are obtained from the first derivative of the transfer matrix [^, |T0|, ^ 



N-l 



H = 2r(,)(E H.,n., + ^irr(o)-i^r'(o) + 



n=l 



trK+{<d) 



(2.27) 
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where K (0) ^ in the second term is introduced for the case when K (0) 7^ 1, r(?7) = snrj 
and the two-site Hamiltonian if„ is given by 

1 

Hn,n+1 = —r^Pn,n+lRn,n+li^) (2.28) 

r[ri) 

Pij is the permutation operator acting on Vi ^ Vj with Vi = Vj = V, i.e. P{x ^ y) = 
y ® X {x,y eV). If 

trir+(0) = 0, (2.29) 



then as pointed out in ref. [0, we will not have a well defined Hamiltonian (2.27) from 



the first derivative of the transfer matrix. In this case, we may obtain a Hamiltonian from 
the second derivative of the transfer matrix. Only when the condition 

tr^[K+{Q)HNom « 1 (2-30) 

holds for the i^-matrix, will such a Hamiltonian be that with nearest-neighbor interac- 
tions. Otherwise, it will contain terms that couple every pair of sites in the bulk with the 
boundary, i.e. the Hamiltonian is non-local |]T0|, |T^. The condition (2.29) holds for the 
/^-matrices given in equations (2.17), (2.18), (2.19), (2.22) and (2.25), but the condition 
(2.30) does not hold for these i^'-matrices. Due to the above reason, therefore, we will 
not discuss the Hamiltonians corresponding to these i^-matrices below. 

Now we focus on the case of trK^{Q) ^ 0. Using the equations (2.2) and (2.8), the 
first term in the equation (2.27), which describes the bulk properties of the XY Z model, 
is given as follows, 

N-l 

H,uik = E [(1 + + (1 - r)a^a^+i + A«^_,], (2.31) 



n=l 



where 



r = ksn 1], A = cnrjdnr]. (2.32) 



and a^,a^,a^ are the Pauli matrices. The second and the third terms in the equations 
(2.27) represent the left and right boundary terms, respectively. Using the equations (2.2) 
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and (2.28), the boundary terms for the ^--matrices (2.20), (2.21), (2.23), (2.24), (2.26) 
and their corresponding "'"-matrices are of the following form. 



boundary 



snr){A_al + B.af + C_a^ + A+a% + B+aj^ + C+a^) 



with coefficients A±,5±,C± being as follows, respectively, 

_ c<-|-dn$-|- 



1± 

B 



± 



C± = 0; 



> {for subcase {h2)) 



B^ = £2"+, B_ ^ a_, 
C_|_ = q;_|_, C_ = £20^-; 



> {for subcase (cl)) 



cn^± dnl;± 



C4 



£20+ 

sn^+ ' 



5_ 



> (/or subcase (c2)) 



5± = t2/x±(A± + 1), C± = T2/i±(A±-l); 



(/or subcase (c3)) 



/I _ cng±d<± 
■ri-± — — . 



2m±(A±-1) . 



> (/or subcase (c4)) 



(2.33) 



(2.34) 



(2.35) 



(2.36) 



(2.37) 



(2.38) 



It is easy to see that the equation (2.20) is the special case = of the equation (2.26). 
Comparing the equation (2.34) with the equation (2.38), we know that the off-diagonal 
elements of the X-matrix are responsible for the B±, C± terms in the Hhoundary 



From the symmetry (A) and the equation (2.27), we know that if the scaling factor f{u) 
satisfies the conditions: /'(O) = and /(O) = const., then f{u) will not change the bound- 
ary terms. f{u) — cnudnu/ (1 — /c^sn^^sn^ii) is such an example. When we take ^_ = /C(/C 
is the half-period magnitude of the Jacobi elliptic function) in the equation (2.23), then 
the equation (2.23) reduces to the equation (2.21) up to a factor cnudnu/{l — fc^sn^w); 
Similarly, the equation (2.26) reduces to the equation (2.24) up to a factor 1/(1 — k'^sn^u). 
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These two factors satisfies the conditions for the scahng factor f{u) above. The equa- 
tions (2.21) and (2.23), therefore, correspond to the same boundary terms upon taking 
= }C,C,^ = — /C; the same holds for the equations (2.24) and (2.26). 



In closing this section, we make a remark on the symmetry of the Hamiltonian (2.27). 
It follows from the equations (2.31) and (2.33) that the Hamiltonian is transposition 
invariant if we make any replacement given below, 

(1) B_^C^, B^^C+- (2.39) 

(2) a)an — ^o-jv-n+i, in particular, ai — > ctat, 0-2 — ^ (yN-i] (2 Af^)) 

b)A^ — > A+, 5_ — > C+, C_ — > B+. 
In the next section, we will find that the second replacement is very useful for the con- 
struction of the Lax pair of the open XY Z spin-chain. 



3 The Lax pair for the open XY Z spin-chain 

In section 2, we have derived the Hamiltonian of the open XYZ spin-| chain. For the ID 
system with sites, the general form of the the Hamiltonian is given by the equations 
(2.27), (2.31) and (2.33)in[, 

N-l 
^ n=l 

+snri{A^al + B^a^ + + A+cx'^ + B+aj^ + C+a^v), (3.1) 

where 

Jx = —2(1 + ksn^T]), Jy = —2(1 — ksn^r]), = — 2cn?7dn?7, (3.2) 
The equations of motion for the system (3.1) are given as follows, 

= -J,<«+l+^n-l) + '^^^n«+l + <-l), (3.3) 

j^al = -J.<«+i + <-i) + </X«+i + <-i), (3.4) 
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|< = + + VnK+i + ^Li), (3.5) 

(n = 2,--.,7V-l) 



d 

Jt 
d_ 

~dt 

d ^ 

df' ^ 



JyGlal + J,a\al + snr7(-2>l_(7? + i{B_ - C_)(7f ), (3.6) 
al = -JXf^2' + '>^x(TX + snr7(2^_(7j'- (5_ + C_)(ji"), (3.7) 
Ja^afa^ + Vi^I - 2isnri{B_at - Caf ), (3.8) 



,v + Jz(Jn-i(jI + sn77(-2yl+< + z(5+ - C+)a^), (3.9) 

^cr^ = -J,<_ia^ + J,,cT^_iCT^ + snr?(2A+a^- (5+ + C+)cT^), (3.10) 



-JxCTN-if^N + Jy'^N-i'^N - '2isnrj{B+ajf - C+a^). (3.11) 



In order to rewrite equations above in the Lax form, we consider the following operator 
version of an auxiliary linear problem, 







1,2,- 


..,N) 


(3.12) 


d T 


= M„(ii)*„, (n = 


2,-- 


.,N) 


(3.13) 










(3.14) 


dt''^ 








(3.15) 



d_ 

di 



where u is the spectral parameter which does not depend on the time t, Ln{u), M„(m), Qi{u) 
and Qn{u) are called the Lax pair. Qi{u) and Qn{u) are responsible for the boundary 
conditions. The consistency conditions for the equations (3.12)-(3.15) are the following 
Lax equations: 

j^Ln{u) = Mn+i{u)Ln{u)-Ln{u)Mn{u), (n = 2, • • • , - 1) (3.16) 

^Lm{u) = Qn{u)Ln{u) - Ln{u)Mn{u), (3.17) 
dt 

^Uu) = M2iu)L,iu)-L,{u)Q,iu). (3.18) 
dt 

Prom the eqs. (3.3)-(3.5) and (3.16), we see that they are the same as those in the case 
of periodic boundary condition. For the XYZ model with periodic boundary condition, 
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the Lax pair L„, M„ has been given by Sogo and Wadati[|T5[. The operators L„ and M„ 
are of the following forms, 



Lr. 



W4 - w^a^ 



(3.19) 



(3.20) 



where 



(3.21) 



(3.22) 



Wi{u) + W2{u) = sni], Wi{u) — W2{u) = ksnrjsnusrL^u + rj), 

Wsi^u) + W4 = sn(M + T]), W4^{u) — W3 = snw, 
and a" = 1, a^,a'^,a'^ = a^^a^^a^ are the Pauli matrices. For the explicit expressions of 
Fj, Gj, Hi see ref. [|^. So, for the open XYZ model, we only have to construct operators 
Qiiu) and Qn{u). If we write Qiiu) and Qn{u) in the following forms: 

Qi{u) = (^7^)2x2, Qn{u) = (aS^^)5 



12x2 



and 



,(-) 



E 



.(+) 



1) 



E4 



Ml I 

O n 



(3.23) 



(3.24) 



1=0 1=0 
Substituting the equations (3.23) and (3.24) into the equations (3.17) and (3.18), we ob- 
tain 32 equations for a-^^'(i, j = 1, 2, / = 0, 1, 2, 3). Due to the following lemma, we only 



have to solve 16 equations ,for example, for a|^'"(2, j 



1,2,/ = 0,1, 2, 3). 



Lemma: If we make the replacement (2.40), then we have 



,(+) 



(3.25) 



provided that the following correspondences hold 



— Arr^ ■ 



(3.26) 
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(3.27) 



(M|)* M3 ; 

under the transposition. Here the superscript t denotes the transposion of a matrix. 

Note that (cr^)* = cr^, (c^)* = — a^, {cr^Y = cr^, then from the equations (3.6)-(3.11) and 
the expressions for M^, M2{i — 1, 2, 3), it is straightforward to show that the conditions 
(3.26) and (3.27) hold. So if we get the solutions for a\'^\ using (3.25) the expressions 
for alj ^ are also obtained. After a lengthy calculations, the expressions for a\f''\i,j — 
1, 2, Z = 0, 1, 2, 3) are as follows, 

a\^' = a{u)b{uy, (3.28) 

= '-^{±{B^-C^)c{u){-2e{u) + f{u)+g{u)) 

-{B± + C^){b'{u) + c'{u) - f{u))}; (3.29) 

+(S± + C±)c{u){f{u) - g{u))}; (3.30) 

(±)3 ^ isnriA± 
ail = G3 — — d{u); (3.31) 

4?° = ^c{u){{B^-Ci)a{u) + {B^ + Ci)b{u)}; (3.32) 

aS?^ = G,±'^^^^b{u){-2e{u)+f{u)+g{u)y, (3.33) 

4^^' = -iG'2±^^^a(x.)(-/(^)+^(^)); (3.34) 

«S?' = T^(/(«)+^(«)){(i?±-C±)a(«) + (5± + C±)6(«)}; (3.35) 

= ^c(^){-(S±-C±)a(^) + (S± + C±)6(i.)}; (3.36) 

4?' = GiT^^^^6(«)(-2e(«) + /(«) + ^(«)); (3.37) 

4?' = iG,±'^^^a{u){-f{u)+g{u))- (3.38) 

4?' = ±^(/(i.)+^(i.)){-(S±-C±)a(^) + (S± + C±)6(i.)}; (3.39) 

4?° = -4?^ (3.40) 
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where 



4?' = '-^{±{B^-C^)c{u){2e{u)-f{u)-g{u)) 

-{B^ + C^){h\u) + c\u) - f\u))}- (3.41) 
4?' = ^-^{{B^-C^){a\u) + c\u)-e\u)) 

T(fi± + C±)c(w)(/H - (3.42) 
^22 = -G?, — — d{u) (3.43) 



a{u) = W2W3 + W1W4, b{u) = W1W3 + W2W4, 
c{u) = W1W2 + w^w^, 

d{u) = {w2-Wif{{w3-Wiy-{w2 + WiY} + 

+ {W2 + wiYiiws + w^f - {W2 - w,Y}, (3.44) 

/ \ 2 2 r/ \ 2 2 / \ 2 2 

The equations (3.28)-(3.43) combined with the operators L„, M„ are the Lax pair for the 
open XYZ spin-chain. In the trigonometric hmit k ^ 0, where sn-u sin-u and taking 
Jx = Jy = ^, Jz = COST], B± = C± = 0, we recover the result in the case of open XX Z 
model given in the ref. [|l3](up to a replacement of 1] by 21]). 



4 Remarks and discussions 

In this paper, based on the analysis of the symmetries of the reflection equations, we 
present their solutions for the open XYZ model case by case. We realised that both 
diagonal and off-diagonal elements of the i^-matrices have contributions to the boundary 
terms of the Hamiltonian, and a specific scaling factor does not affect the properties of 
the system. Furthermore, we construct the Lax pair for the open XYZ model to show 
its integrability in another way. We believe that our procedure to the solutions of the 
reflection equations can be applied to other problems. Recently, the classification of six- 
vertex, eight- vertex solutions of the coloured Yang-Baxter equation has been given [IB, pT[ , 
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the Yang-Baxter equation with dynamical parameters (or called the Gervais-Neveu-Felder 
equation) has been studied[|l8[, we expect to investegate their corresponding reflection 
equations on the basis of analysing their symmetries. 
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